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A numerical program to calculate the photoproduction of prompt photons is presented. The 
code includes full next-to-leading order corrections to all subprocesses. The results are com- 
pared to recent ZEUS data. 



1 Introduction 

High energy electron-proton scattering is dominated by photoproduction processes, where the 
electron acts as a source of quasi-real photons which interact with the partons in the proton. 
These processes are of special interest since they are sensitive to both the partonic structure of 
the proton as well as of the photon. In particular, they offer the possibility to constrain the 
gluon distribution in the photon, since the subprocess qg — > jq, where the gluon is stemming 
from a resolved photon, is contributing already at leading order. 

The first calculationscla of higher order corrections to the Compton process ^q — ► ^yq con- 
sider only the fully inclusive cross section without the possibility to deal with isolated photons. 
More recent calculations done by Gordon/VogelsangtJ and Krawczyk/ZembrzuskiH take isolation 
into account, but only by adding a subtraction term evaluated in the collinear approximation 
to the fully inclusive cross section. Moreover, the code ofH does do not contain the full set of 
NLO corrections. 

The calculation presented here takes into account the full NLO corrections to direct as well 
as resolved and fragmentation parts. In addition, it includes the box contribution g'j — > gj 
which is formally an 0{a 2 s ) correction, but known to be importantH. A major advantage of the 
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present code is also given by the fact that it is constructed as a "partonic event generator" and 
as such is very flexible. It produces N-tuples of partonic final state configurations which can be 
generated once and for all. Based on these N-tuples, suitable observables matching a particular 
experiment can be defined and histogrammed. 

2 Theoretical formalism 

The inclusive cross section for ep — ► 7X can symbolically be written as a convolution of the 
parton densities of the incident particles (resp. fragmentation function for an outgoing parton 
fragmenting into a photon) with the partonic cross section a 

da ep ^ x (P p ,P e ,P y ) = J2 I dx ej dx pj dzF a/e (x e ,M)F b/p (x p ,M p ) 



a,b,c 

da ab ^ cX (x p P p , x e P e ,P^/z, fx, M, M p ,M F ) D l/c {z, M F ) (1) 

where M, M p are the initial state factorization scales, Mp the final state factorization scale and \x 
the renormalization scale. The subprocesses contributing to the partonic reaction ab — > cX can 
be divided into four categories which will be denoted by 1. direct direct 2. direct fragmentation 
3. resolved direct 4. resolved fragmentation. The cases 1. and 3. correspond to c = 7 and 
D„/ c (z,Mf) = 5 n 5(l — z) in ([!]), that is, the prompt]] photon is produced directly in the hard 
subprocess. The cases 1. and 2. correspond to a = 7, with F 1 / e approximated by the Weizsacker- 
Williams formula for the spectrum of the quasi-real photons 
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For the resolved contributions (3. and 4.) a parton stemming from the photon instead of the 
photon itself participates in the hard subprocess, such that F a / e (x e , M) is given by a convolution 
of the Weizsacker- Williams spectrum with the parton distributions in the photon: 

F a /e{xe,M) = J dydx-y f^(y)F a / J (x 1 ,M)5(x 1 y - x e ) . (3) 

At next-to-leading order, the 0{a s ) corrections to the corresponding subprocesses are taken 
into account. The initial state collinear singularities are absorbed into the functions i ? a / 7 (x 7 , M) 
resp. F b /p(x p , M p ) at the factorization scales M,M p , the final state singularities are absorbed 
into the fragmentation functions D^/ C (z, Mp) at the fragmentation scale Mp. As a consequence, 
the four subparts separately depend strongly on M, M p , Mp and only the sum of all four parts, 
where the leading scale dependence cancels, has a physical meaning. 

From a technical point of view, there are basically two methods to isolate the infrared 
singularities appearing in the calculation at NLO: The phase space slicing method and the 
subtraction method. The method used here is basically a phase space slicing method. For 
further details see El. 

In order to single out the prompt photon events from the huge background of secondary 
photons produced by the decays of 7r°,r/, u mesons, isolation cuts have to be imposed on the 
photon signals in the experiment. A commonly used isolation criterion is the following: A 
photon is isolated if, inside a cone centered around the photon direction in the rapidity and 
azimuthal angle plane, the amount of deposited hadronic transverse energy El T ad is smaller than 
some value E™ ax fixed by the experiment: 

(V ~ V,) 2 + (0 " 7 ) 2 < R* , E^ ad < ET X (4) 



By "prompt" we mean that the photon is not produced from the decay of light mesons. 



Following the conventions of the ZEUS collaboration^, we used EVp ax = 0.1 pj, and R exp = 
1. Isolation not only reduces the background from secondary photons, but also substantially 
reduces the fragmentation components. 



3 Numerical results and comparison to ZEUS data 

The numerical results can be presented only briefly here, for further details the reader is referred 
to @. For the parton distributions in the proton the MRST2 parametrizationi is taken. The 
default choice for the parton distributions in the photon is AF(fl, for comparisons we also used 
the GRV0 distributions transformedto the MS scheme. For the fragmentation functions we use 
the parametrization of Bourhis et alO. We take nj = 4 flavors and for a s (/j,) we use an exact 
solution of the two-loop renormalization group equation, and not an expansion in log(/x/A). 
The scales have been set equal to ; a variation of the scales between p^/2 and 2pj, leads to a 
variation of the cross section by less than 10%. The rapidities refer to the ep laboratory frame, 
with the proton moving towards positive rapidity. To match the ZEUS conventions, we set 
Qmax = lGeV in the Weizsacker- Williams spectrum and restrict the photon energy y = E^/E e 
to the range 0.2 < y < 0.9. 




Figure 1: (a) Comparison of NLO to LO result for the photon rapidity distribution, with isolation, (b) Effect of 
changing the lower bound on y. Solid line: default, dotted line: lower bound on y decreased by 10%, dashed line: 
0.2 < y < 0.9 with GRV photon structure functions, (c) and (d): Comparison to 1996/97 ZEUS data. 



Figure ||(a) shows a comparison of the NLO to the leading order result and displays the 
magnitude of the box contribution. The higher order corrections enhance the isolated cross 
section by about 40%. Fig. |l](b) illustrates the sensitivity of the cross section to a variation of 
the energy range of the photon: A 10% change of the lower bound on y leads to a spread which 
- except in the low rapidity region - is larger than the one caused by a different set (GRV) 
of photon structure functions. Note that experimentally, the energy of the incoming photon in 
photoproduction processes is reconstructed from the final hadron energies. In order to obtain 
the "true" photon energy y, corrections for detector effects and energy calibration have to be 
applied. Fig. |l](b) shows that a good control on the error in the reconstruction of y is crucial 
for a detailed comparison between data and theory. 

Figs, [l] (c) and (d) show a comparison to ZEUS 1996/97 dataQ for the photon pr and rapidity 
distributions with two different sets of photon structure functions ( AFG and GRV) . For the pt 
distribution the agreement is quite satisfactory; in the rapidity distribution there is a tendency 
that theory underpredicts the data in the backward region. 



4 Conclusions 



A full NLO program for the photoproduction of prompt photons has been presented. The code 
generates N-tuples of partonic final state configurations which serve as a starting point to define 
appropriate observables matching the experimental situation. The agreement between ZEUS 
1996/97 data and theory is in general satisfactory; a discrepancy can be observed at low photon 
rapidities. With the present errors on the data, a discrimination between the AFG/GRV sets of 
parton distributions in the photon is not possible, but a forthcoming high statistics analysis of 
all 1996-2000 data announced by the ZEUS collaboration will improve this situation. 
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